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Abstract 

A treatment of the modal decomposition of the pressure field in a combustor as deter- 
mined by two Kulite pressure measurements is developed herein. It is applied to a Pratt 
& Whitney PW4098 engine combustor over a range of operating conditions. For modes 
other than the plane wave the new part of the treatment is the assumption that there are 
distinct frequency bands in which the individual modes, including the plane wave mode, 
overlap such that if circumferential mode m and circumferential mode m — 1 are present 
than circumferential mode m — 2 is not. Consequently, in the analysis used herein at fre- 
quencies above the first cut-off mode frequency, only pairs of circumferential modes are 
individually present at each frequency. Consequently, this is a restricted modal analysis. 
A new result is that the successful use of the same modal span frequencies over a range of 
operating conditions for this particular engine suggests that the temperature, T, and the 
velocity, v, of the flow at each operating condition are related by c 2 — v 2 = a constant where 
c is the speed of sound. 


Nomenclature 


A amplitude of clockwise pressure wave 
B amplitude of counter-clockwise pressure wave 
B e resolution bandwidth, Hz., B e = 1/Tg = r/NP = 11.71875-ffz 

D delay time r, sec. sometimes expressed as a number of samples, D = t 48000) 
/ frequency, Hz 

f c upper frequency limit, f c = 1/2 At = r/2, Hz. (24000 Hz.) 

G xx (f ) power spectral density 
G xy (f) cross power spectral density 
j positive imaginary square root of — l,v/— T 

L c cost function 

L y number of frequencies, f c /Af = A/2 ( 2048 ) 

M Mach number, M = uc 

m mode number running index 

n mode number running index 

n 0 e number of overlapped segments, with 50 % overlap n 0 e = 468 

NP segment length, number of data points per segment (4096) 

rid number of disjoint (independent) data segments/blocks,? - ^ = B e T to t a i = 234 

no contaminating noise at station 0 

ng contaminating noise at station 6 

P(f) Fourier transform of p(t) 

p(t) pressure signal 

Pk conversion constant, lO*- 171 / 20 ^ /0.25 
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pit) pressure signal 

r sample rate, samples/sec. (48000) 

i?(r) radial component of the pressure 

51 signal at station 1 

5 2 signal at station 2 

Ttotai total record length, sec. (« 20 sec.) 

t time 

Td(i) record length of segment i 

u velocity vector 

Subscripts 

0 microphone reference location 

1 running segment index 

Im Imaginary part 

Re Real part 

x signal x{t) 

y signal y(t) 

Symbols 

A / frquency step, l/T^Hz. (11.718) 

At sampling interval, 1/r (1/48000), sec. 

7 2 n analytical coherence threshold 

7 2 y coherence function 

7 (/) estimated magnitude squared coherence (MSC) function 
u) weighting function also known as a window function 
c t> phase angle 

9 angular microphone location 

Ti propagation time delay 

Superscripts 

* complex conjugate 

+ clockwise rotating mode 

— counter-clockwise rotating mode 

ensemble average 


I. INTRODUCTION 

To understand combustion noise measured in the far field of the PW4098 engine it is necessary to 
investigate and understand the pressure field in the annular combustor. The data analyzed is from a Pratt 
& Whitney PW4098 engine. The test was conducted as part of phase 2 of the NASA Engine Validation of 
Noise Reduction Concepts (EVNRC) Program. For this test two Kulite pressure transducers were mounted 
in the combustor and four far field microphone were used. The coherent combustion noise seemed to have a 
modal pattern. Consequently, a restricted modal model was developed and tested to see if it was compatible 
with the available data. This work is reported herein. To be useful, the pressure field must be measured and 
analyzed in a manner consistent with the physics of the propagation process. A scheme that accomplishes 
this is discussed herein. 

The Kulites failed during the test. Consequently, a sensor validation analysis was conducted by Miles. 1 
Results of this study were used to select the data evaluated herein. The Kulite signals and the far field 
microphone signals are used to obtain the far field coherent combustion noise output power and estimate the 
amount of coherent combustion noise that appears in the far field. 2 A method for separating correlated noise 
sources (core noise) and uncorrelated noise sources (jet noise) in far field measurements of turbojet engine 
core noise using multiple microphones has also been developed and tested 3 .Results from the Kulite validation 
study 1 and this study of annular duct modes measured in the combustor by Kulite pressure transducers is 
used to interpret results of Miles. 2,3 

The pressure field inside an annular combustor is governed by the solution of an eigenvalue equation. 
The pressure is assumed to propagate as if it were in an infinite duct. Consequently, boundary conditions 
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at the duct ends, reflections from the duct ends, and the possibility of standing waves are not considered. 
The pressure does propagates in particular bands identified by frequency ranges. The pressure pattern in 
a combustor is analyzed by breaking it up into the natural frequency bands of pressure propagation in 
an annular duct. Frequencies in the lowest band propagation in a plane wave mode. Higher bands move 
in propagating modes composed of eigenmode solutions to the annular duct eigenvalue equation governing 
wave propagation in an infinite annular duct with wall boundary conditions. The modes are composed 
of sinusoids in a circumferential direction and combinations of Bessel functions in a radial direction. The 
physics of pressure waves propagating in annular ducts is discussed by Tyler and Sofrin 4 in a treatment 
of axial flow compressor noise and is presented by Morse 5 as a problem on page 603. Each higher mode 
propagates only above its own cut-off frequency. It is expected that the pressure field in a combustion duct 
extends over a frequency range in which several modes are propagating. The distribution of energy in the 
frequency band for each propagating mode is unknown and might provide insight into the physics of the 
propagation process. Consequently, our understanding of the pressure field in an annular combustion would 
be increased if one could determine the distribution of energy between duct modes. 

The use of normal mode theory to analyze the propagation of noise in ducts to find the optimal microphone 
positions in order to estimate sound power from sound-pressure measurements was studied by Dyer. 6 This 
analysis was soon extended to look at spinning acoustic modes generated by fans. 4, 7-11 This analysis has 
also been used to look at the modal content of noise generated by a jet in a pipe. 12-15 An acoustic modal 
analysis of a YF102 combustor installed in a ducted test rig was conducted by Karchmer 16 who used the 
assumption of equal amplitude clockwise and counter-clockwise spinning modes in his analysis. He found 
that he only needed the circumferential modes 1 through 6 and the corresponding zeroth order radial modes 
in his reconstruction of the measured data. Acoustic modal analysis of the pressure field in the tailpipe 
of a turbofan engine was conducted by Krejsa and Karchmer 1 ' again the assumption of equal amplitude 
clockwise and counter-clockwise spinning modes was used in the data analysis. 

Before duct mode theory was applied to study core noise other techniques were used. Core noise from 
a Pratt & Whitney JT8D was studied by Grande 18 using multiple microphone measurements to obtain 
auto and cross power spectra in a tail pipe extension as well as far field auto power spectra. Ten internal 
microphones were used to obtain measurements of narrow-band and one-third octave band pressure level 
spectra in a Avco Lycoming YF102 combustor by Wilson and 0’Connel. 19-21 Core noise in General Electric 
engines has been studied by Matta, Sandusky, and Doyle 22 to relate performance with emissions and noise. 
Coherence functions and transfer functions were used by Doyle and Moore 23 in a core noise investigation 
of the CF6-50 Turbofan Engine. They used five Kulite pressure transducers in the combustor, nine other 
internal sensors, and 15 far-field microphones. Doyle and Moore 23 removed the time delays between the 
internal and far field signals using cross-correlation analysis. Doyle and Moore 23 set the number of averages 
they used to 100 and as a consequence data below coherence function values of less than 0.1 are ignored. In 
this paper, the coherence function is plotted on a log amplitude scale since the coherence floor is near 0.01. 
Doyle and Moore 23 try to identify the source location from vectoring of cross-correlation time delays and 
do not do a model study. However, examination of internal coherence and transfer functions plots indicate 
that combustion modes might be present. In addition, Doyle and Moore 23 do state that the appearance of a 
double peak with positive and negative time delays in some of the cross-correlation plots suggest that waves 
are moving circumferentially around the combustor. 

The problem of mode propagation in an annular combustor is complex due to the presence of high speed 
flows and temperature gradients. The results given here are limited to using circumferential modes which 
resemble plane waves and the plane wave mode to replicate the measurements. To analyze the pressure 
distribution in a duct into its modes, it is necessary to be able to measure the relationships among the 
pressures at many different points. However, for this test measurements were available at only two points. 
Consequently, in this paper we only investigate the feasibility of using a restricted pressure modal model to 
replicate the single measured relationship. 

The investigation discussed will: 

• Show the Kulite instrumentation was functioning. 

• Provide evidence the peaks and dips in the measured auto-spectra and cross-spectra magnitudes are 
due to duct propagation modes rather than a fluctuating pressure generating mechanism. 

• Identify the energy distribution in the plane wave mode. 
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PW4098 Engine at C-ll Stand for EVNRC Phase 2 Tests 



Figure 1. Pratt &: Whitney test stand Cll, West Palm Beach Florida, with PW4098 engine and attached 
acoustic inflow control device also with and without aft acoustic barrier walls for EVNRC Phase 2 tests. 
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• Examine the validity of the assumption that one has equal amplitude clockwise and counter-clockwise 
spinning modes. 

• Test the idea that the data can be reconstructed using the first few circumferential modes, m = 1, m = 
2, m = 3... and the corresponding zeroth order radial modes p = 0. 

• Determine the energy distribution in each band of frequencies composing a propagation mode 

II. Experiment 

The two Kulite pressure measurements made in a 
Pratt & Whitney PW4098 combustor will be dis- 
cussed next. The measurements were made in a 
study of aircraft engine core noise conducted as part 
of the NASA Engine Validation of Noise Reduction 
Concepts (EVNRC) Program. The Kulite at 127° is 
identified as Kulite 1 and produces a pressure signal, 

Po(t) while the other at 337° is identified as Kulite 
2 and produces a pressure signal, pe(t). Kulite an- 
gles are measured clockwise from top dead center 
viewed from the rear. The combustion chamber is 
annular. The PW4000 was fitted with 24 injectors. 24 
Fig. 1 on the preceding page shows the test stand. 

The spectral estimate parameters are shown in ta- 
ble 1 on the next page The signal processing algo- 
rithms used were written in Fortran. They are based 
on subprogram modules developed by Stearns and 
David. 25 They were modified for this project to pro- 
vide for time delay selection. In the calculations 
the segments were overlapped by 50 percent. Fig. 2 
shows the Kulites mounted in the combustor with the water cooled jackets. Ten sets of data were analyzed. 
This paper shows results from three sets of combustor measurements made on the first test day. The values 
of N1 CORR. used are: 1622 rpm, 1750 rpm, and 1900 rpm. 

III. Analysis 

A. Measured Auto-Spectra and Cross-Spectra 

For discrete time signals the computation procedure is as follows. The total record length T to tai for received 
signals pe(t) and po(t) is divided into rid disjoint (independent) data segments/blocks, each of length Td = 
Ttotai /nd so that there is an ensemble of measurements {pei(t)} and {poi(i)} 



Figure 2. Kulite hardware mounted on PW4098 tur- 
bofan engine at Pratt &; Whitney test stand Cll, West 
Palm Beach Florida for EVNRC Phase 2 tests. 


{P 9 i{t ) } = {s 2 (t) + n e }-, * = 1, • • • ,n d 0 < t < T (1) 

{rhM} = {si(t) + noWlj , i = l,---,n d 0 <t<T 

where ng(t) and no(t) are the contaminating noise at receiver location 1. Angular positions, 8 , are relative 
to Kulite 1 at 227° which was used as an arbitrary reference 0 = 0° location. The relative angular location 
of Kulite 2 at 337° is then 9 = 210° . 

Using notation similar to that used by Piersol 26 the one sided auto-spectra, Goo(f),Ggg(f) and the cross- 
spectrum, Goe(f) are estimated for a given pair of random pressure signals, at a frequency f as 


Goo(f) 


Po(f)P 0 *(f) 

2 nd 

5^EW)Po.(/1 

2=1 


( 2 ) 
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Gee(f) 


Goe(f) 


Po(f)PSif) 

o n d 


T d n d 


E P e*C f)Pei(f) 


Po(f)P 9 *(f) 

2 


T d n d 


T, p ei(f) p 0i(f) 


i—1 


( 3 ) 

( 4 ) 


where the “ indicates an ensemble aver- 
age and Pgiif) and Poi(f) are the Fourier 
transform of the measured pressure signals 
at the reference position and at the angle 
9. The * designates the complex conjugate. 
The Fourier transforms used are windowed 
Fourier transforms: 



where ui{tau) is a weighting function also 
known as a window function. This approach 
is used in evaluating spectra of random data 
(see Stearns 27 section 15.6). In practice 
these quantities are computed using a Fast 
Fourier transform. The estimate Goe(f) is a 
complex number such that 


Table 1. Spectral Estimate Parameters 


Parameter 

value 

Segment length i.e. Data points per segment, NP 

4096 

Sample rate, r, samples/second 

48,000 

Segment length, T d = NP/r, seconds 

0.08533 

Sampling interval, At = 1/r , seconds 

2.0833 x 10“ 5 

Frequency step, A / = 1/Td, Hz 

11.718 

Upper frequency limit, f c = 1/2A t = r/2, Hz 

24000 

number of frequencies, Ly = / C /A/ = NP / 2 

2048 

Time delay, r = 6323/48000, seconds 

0.1317 

number of independent samples 

234 

overlap 

0.50 

Sample length, sec. 

20 


Goe(f) = G Re (f) - =| G 0B (f) \ e ~^ .(/) 

where 


( 6 ) 


I GM I 

0o e(f) 


—2 


—2 


G Re (f) + G Im (f) 
-i G Im {f) 


1/2 


tan 


[G Re (f)\ 


( 7 ) 


From Piersol 26 and Stearns and David 25 the spectral estimate will be at discrete frequencies separated by 
A / = 1/T d . For digital data processing with a sampling rate, r, and NP data points per segment, the length 
of each segment is T d = N/r seconds. The sampling interval is At = 1 jr. Consequently, A / = 1/(A tNP) 
and the upper frequency limit is f c = 1/(2 At) = r/2. The, number of frequencies for nonredundant spectral 
estimates is L y = / C /A/ = Nj 2. The spectral estimate parameters are shown in Table 1. 

The signal processing algorithms used were written in Fortran based on subprogram modules 25 which were 
modified for this project to be in Fortran 90 form. In addition, the processing scheme was changed to permit 
delaying one time sequence with respect to another. For these computations, a data segment/block record 
length of T d = 4096/48000 = 0.08533sec. seconds was used and the number of disjoint data segments /blocks 
was n d = 234. Consequently, total record length was T d n d ~ 20 sec.. The analysis was done with 50 per 
cent over lap so n 0 i ~ 468. 

The coherence is given by 


7oe(/) 


\Goe(f)\ 2 

G 00 (f)G ee (f) 


( 8 ) 
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The measured coherence calculated using segments overlapped by 50 % is given by 


—2 m IE^(/W/)I 2 M 

lo0[J> Er=i \ p oi(f )\ 2 Er=i \ p oi(f )\ 2 u 

Given measurements from two sensors to examine, we can say a sensor is providing an incorrect signal if 
the signals are uncorrelated with each other or the correlation changes abruptly in some irregular fashion. A 
procedure developed by Miles 1 based on the concept of aligned and unaligned coherence is used to determine 
the validity of the signals from the Kulites. Miles 1 shows that by comparing a coherence function calculated 
using aligned and unaligned time histories one may decide if the signals from two sensors are uncorrelated 
or have changed in some irregular fashion. The Kulite pressure data set examined in this paper was selected 
using this procedure. This procedure compares the normal measured coherence, 7 q g(f), called the aligned 
coherence, with one calculated by time delaying one signal by an amount D greater than the segment /block 
interval Td . The value of D is chosen so that the two signals are not in the same segment/block interval, 
Td and consequently appear to the processing procedure to be uncorrelated. For this study a value of D 
approximately 1.5 times the segment/block interval is used. The measured unaligned coherence, 7 q e (/, D) is 
not zero. Even if no tones are present, the coherence of two disjoint (incoherent/independent) random noise 
records has a value dependent on the number of independent data segments/blocks used to calculate the 
coherence, rid ■ In Miles 2 results were obtained from computer simulation that show good agreement with 
the theoretical estimate of the analytical coherence threshold 


Inn = 1 — (1 — ( 10 ) 

where we use P = 0.95 and rid = n 0 e. The coherence threshold y 2 n is discussed by Carter, 28,29 Halliday 
et. al. 30 (page 247), and Brillinger 31 (page 317). The coherence threshold y 2 n has a value which is 
greater than 95 % of the values of the coherence of two independent time series calculated using rid disjoint 
data segments/blocks. For the values shown in table 1 on the preceding page with a 50% overlap, a total 
record length, TTotai, of 20 seconds, and a resolution bandwidth, B e , of 11.718 Hz., the number of data 
segments/blocks is rid = ZAfTTotai ~ 468. The coherence of random noise for this case is 7 2 n = 0.00639431. 
The concept of aligned and unaligned coherence is discussed in more detail in Miles. 2 

In addition, if tones are present in the time signal they will be present in all ensembles averaged and 
appear in the unaligned and aligned coherence. In addition to aiding the identification of persistent tones, the 
unaligned coherence provides a reference coherence. Small values of the aligned coherence greater than the 
unaligned reference coherence can be easily identified as significant since one knows what the uncorrelated 
coherence looks like. However, the major support it provides is in determining if a sensor is providing a 
useful signal. If the aligned and unaligned coherence are similar than the two sensors are uncorrelated and 
one sensor has failed. 

For the three cases discussed herein measured auto-spectra magnitudes and cross-spectra magnitudes 
and phase angles are shown in Fig. 3 on the next page . Both a measured and an unwrapped phase 
angle are shown. While the amplitudes of the cross-spectra are not remarkable, the cross spectral phases 
are quite notable in that the phase differences are not as predominantly either 0° or 180° as was found by 
Karchmer 16 in his study of YF102 combustor modes. The presence of this feature in the YF102 phase spectra 
which is a characteristic of standing circumferential waves was found at all conditions tested by Karchmer. 
Consequently, in the PW4098 combustor the circumferential waves are not predominantly standing waves 
but are spiraling or spinning waves. The amplitudes of the cross-spectra and the coherence function do have 
a lumpy appearance which is characteristic of the presence of a modal structure. 

In Bendat 32 and in Piersol 26 the random error in the phase estimates due to statistical sampling is given 
in terms of the standard deviation of the estimated phase angle, 9 oe, by 


o- [0o e (/)] ~ sin 1 


[i-i,y/)] 1/2 l 

|7oe|V2rw J 


( 11 ) 


where a [0 o e(/)] is measured in radians and as used herein n 0 g is selected to be the number of overlapped 
segments or blocks used in the spectral calculations. For the special case where the term in curly brackets 
is small Eq. 11 becomes 
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Figure 3. Auto-spectrum magnitude for Kulite (1) at 127 degrees and Kulite (2) at 337 degrees and cross- 
spectrum magnitude and phase between Kulite (1) and Kulite (2). 
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l7oe| V2n 0 £ 


( 12 ) 


where for the unknown coherence 7 q S (/) the estimated coherence 7oe(/) from Eq. 9 is used. A plot of the 
standard deviation of the phase angle in degrees verses coherence is shown in 4 for n 0 t = 468. 

When the coherence is greater than 0.02, Fig. 3 shows the standard deviation should be less than 13 
degrees for n 0 i = 468. This is accurate enough if no spinning modes are present and the phase angle difference 
is either 0° or 180°. 



Figure 4. Standard deviation of phase angle of Gqq based on 7 ^ and n 0 < . 
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The measured aligned and unaligned coherence is shown in Fig. 5 for the three cases. While many 
coherence values are less than 0.1 the coherence is still valid. As discussed, the aligned coherence is grater than 
the unaligned coherence indicating the coherence is not due to random noise. The random noise coherence 
value is the unaligned coherence value. In addition, the measured unaligned coherence is generally below the 
95% confidence value for the analytic coherence threshold of the random noise which is = 0.00639431 
and is shown as the black line in Fig. 5(a) through 5(c). Fig. 5(c) shows that at this 2600 rpm test point 
near the end of the test run the measured coherence is the same as the random noise coherence. Looking at 
Fig. 3(c) we see that Kulite 2 has failed and the measured auto spectrum G 22 has an incorrect value. Using 
a statistical coherence threshold procedure discussed by Miles 1 each test point was examined for validity and 
sensor degradation. Initially runs with rpm in the range of 1622 to 2400 were used in the modal analysis 
However, eventually runs with 2304 rpm and 2400 rpm were discarded. This left eight cases. 



FREQUENCY, Hz FREQUENCY, Hz 

(a) 1622 rpm ( N1 Corr ) (b) 1999 rpm ( N1 Corr ) 



0 200 400 600 800 

FREQUENCY, Hz 


(c) 2600 rpm ( N1 Corr ) 


Figure 5. Comparison of aligned and unaligned coherence between combustor Kulite (1) at 127 degrees and 
combustor Kulite (2) at 337 degrees. 


B. Modal Decomposition Method 

The approach implemented here is based on the mathematical modal decomposition technique used by 
Karchmer. 16 Technical details on the eigenvalue problem and solution are given in Appendix A. The problem 
solved includes the mean flow. A new result given in Appendix A is that the successful use of the same modal 
span frequencies over a range of operating conditions for this particular engine suggests that the temperature, 
T, and the velocity, v , of the flow at each operating condition are related by c 2 — v 2 = a constant where c 
is the speed of sound. The use of modal span frequencies will be discussed in the next section. 

The analysis starts by assuming that a randomly occurring instantaneous pressure disturbance produced 
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by the combustion process sends a clock-wise and counter clockwise pressure wave moving around the annulus. 
At an angle 8 the two waves add to produce the resulting pressure written in terms of annular duct modes: 


Pe(t) = P+(t) + P e ~(t) (13) 

M—l 

= ejWt Y + B rn eA^ +m9 '> 

m—0 

where A m and B m are the amplitudes of the clock-wise and counter-clockwise pressure waves associated 
with the lowest radial order of the mth circumferential mode, <p m is the phase of the mth mode, j = yj— 1 
and M is the total number of modes present. Note higher order radial modes are neglected. 

At the reference angle where 8 = 0 


Po(t) = P+(t) + P 0 ~(t) (14) 

M—l 

= J2 K* + B ™ } e j * m 

m—0 

Consequently, the cross-spectrum PqPq is 

M—l 

P 0 Ps = Y + B "*] e ^ m 

m= 0 
N—l 

Y \A n e~A^- ne ) + B n e-A^+ ne ) 

n— 0 

M—l N-l 

= Y E + B mB n )eA^-^+ nS ) 

m= 0 n— 0 
M—l N- 1 

+ E ( BmBn + A n B m )eA^-^- nB ) (15) 

m= 0 n— 0 

After taking the ensemble average the products become correlations. The assumption is now made that the 
modes are independent. Consequently, the cross-correlations vanish. All terms in the averaged version of 
Eq. ( 15 ) for which m ^ n vanish. 

M - 1 

p^p! = Y ( A ^ +B ^) ejm0 

m—0 

M—l 

+ E + AA) e ~ jme ( 16 ) 

m—0 

Using the definition e’ m6 = cos (md) + j sin (md) and noting that 

(A m + B m ) 2 = A^+2B^+Bl (17) 

it follows that Eq. (16) becomes 

M—l 

PoPg = Y, (Am + B m ) 2 cos m8 

m—0 

M—l 

+3 Y ( A2 m- B h ) sin md 

m—0 
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B oP; Re +jP 0 p; Im 


\PoPg\e^ oe 

(18) 

[(PoPgAe) + (PoPgflml V2 

(19) 


Equation ( 18 ) is the general expression for the complex cross-spectrum between the pressure measured at 
the reference position and the pressure measured at a given angle 9 relative to the reference position. The 
magnitude is given by Eq. ( 19 ). 

The phase of this cross-spectrum , then, is 


Ml / 1 TD \ 4 .„„— 1 { B oPg )lm 

(PoByAm, B m ) — tan 

(PoPg)Re 

. Em=o ( A ™ ~ B m) sin m 9 

= tan -1 v ^ 

Em =0 ( A ™ + B ™) 2 COSTO0 


( 20 ) 


Table 2. Mode span 


Mode, m 

Low index 

High index 

Low frequency, Hz 

High frequency, Hz 

0 

0 

15 

0 

175.78125 

1 

16 

28 

187.50000 

328.12500 

2 

29 

46 

339.84375 

539.06250 

3 

47 

71 

550.78125 

832.03125 

4 

72 

73 

843.75000 

855.46875 

5 

74 

75 

867.18750 

878.90625 

6 

76 

77 

890.62500 

902.34375 


C. The restricted modal model 

Equation ( 18 ) and ( 20 ), together with the single cross-spectrum measurement from the PW4098 combustor, 
form the basis of the modal analysis feasibility study conducted herein. The approach used by Karchmer 16 
assumed that only the first six modes make a significant contribution to the pressure field in the combustor 
over the frequency range of interest in core noise studies. Consequently, to determine six unknown values of 
A m and six unknown values of B m one would use measurements of cross-spectrum magnitude and phase at 
twelve positions ( 9 ) to create twelve equations so that one could solve for twelve unknowns. 

However, in the study of the YF102 combustion modes Karchmer 16 found that the measured cross- 
spectrum phase was either 0° or 180° depending on the sign of the real part of the cross spectrum. Conse- 
quently, from Eq. ( 20 ) the imaginary part of the cross spectrum is zero and in evaluating modal constants 
A m and B m he could assume A m = B m . Consequently, he could then infer that the sources are such 
that an instantaneous pressure disturbance created in the combustor sends equal amplitude clockwise and 
counter-clockwise pressure waves traveling circumferentially around the combustor annulus. This creates a 
standing circumferential pressure wave rather than a rotating or spinning circumferential pressure wave. For 
this case, the Kulite cross-spectra phase angles would have only values of 0 degrees or 180 degrees. Using 
this information one only needs measurements of cross spectrum magnitude at six position (9) to create six 
equations so that one can solve for six unknown values of A m where A m = B m . Consequently, the first 
assumption made herein is that for this case the pressure is a standing wave and A m = B m . This is done 
even though some of the phase angles measured seem to have intermediate values. 

The theory of mode propagation and the measurements by Kerschen and Johnston indicate that modes 
cut-on at particular frequencies. 12-15 Consequently, the second assumption made herein is that one can 
choose cut-on frequencies for each mode used to model the system. The selected mode span frequencies are 
shown in Table 2. 
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Examination of the modal decomposition done by Karchmer 16 shown in Fig. 7 in his paper indicates that 
one may easily assume that only the most recent cut-on mode is important. Consequently, in this model when 
mode m is cut-on then the amplitude of the m — 2 mode is set to zero. At any frequency only the amplitude 
of mode m and mode m — 1 must be found. This will be designated using the notation m _ 1 . Due to this 
restriction on the modes, the term restricted modal model is used herein. The final assumption is that the 
reconstructed model auto spectrum, |P 0 PJ| 2 should be less than or equal to the measured auto-spectrum. 

Using these assumptions we have a measurement of cross spectrum magnitude and phase and two equa- 
tions and two unknowns, A m (f) and A m _i(/) at each frequency. The two equations are written in complex 
notation as 


\G 0 e(f)\e j ^ f) = |P 0 P;(A m (/),A m _ 1 (/))|e^^(^’ A ™- 1 ^» (21) 

where 


P 0 P g *(A m (f),A m ^(f)) = 
| P 0 P 9 *( A m (/), A m _ 1 (/))| 2 = 


<f>oe(f) = 


| P 0 P 0 *( A m (/), A m _ i (/))| 2 = 

where a mode m is selected to start at a particular frequency. 

The assumptions made that only mode m and mode m — 1 are active at a frequency, that a correct 
particular frequency can be selected for a mode to cut-on, and that the pressure wave is stationary and not 
spinning are not perfectly valid. In addition, the measured cross-spectrum required as input to the two equa- 
tions is experimentally determined and subject of nominal experimental error and statistical uncertainties. 
Also, the measured auto-spectrum includes random noise in addition to coherent signals from propagating 
waves. Consequently, a solution method that provides an optimum solution in a least squares sense without 
derivatives was used. Algorithms for minimization without derivatives are discussed by Brent. 33 The search 
technique used in this study is described by Powell 34 and Fortran computer code is given in Shapiro 35 and 
Kuester 36 . The code used was a modified version of the one given by Shapiro 35 which was updated to be in 
a FORTRAN 90 style. The cost function used is written in terms of sound pressure level and phase angle. 
We have 


4 A^ cos (m0) 

m,m — 1 

| E 4 ^m COS ( m6 0 

\m,m— 1 y 


4*00 (A-m-, A-m—l) — 


0° or 360° : P 0 P$ > 0 
180° : Pjf < 0 


E 


( 22 ) 

(23) 

(24) 

(25) 


Loe(f) 

— 10Logio(\Goo{f)\ 2 ) 


(26) 

L X x(f) 

= min 10Aogio(|Gn(/)| 2 ), 10Logi 0 (|G 2 2(/)|^) 


(27) 


= 10log 10 (\P 0 Pg (A m (/), A m _i)| 2 ) + 20log w {p k ) 


(28) 

Hoo(m,f) 

— 10Zogio(|Po-Po*(A m (/), A m _i)| 2 ) + 20logio(pk) 


(29) 

0o e(m,f) 

= 0oe(A m (f), A m _i(/)) 


(30) 

Cl(f) 

= [Loe(f) ~ H og (mJ)] 2 


(31) 

G2(/) 

f 0 : H 00 (m,f)<L xx 

1 [Lxx(f) - H 00 (mJ)} 2 : H 00 (m, f) > L xx 

(/) 

(/) 

(32) 
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C3(/) 


C 


min [(0 oe (/) - 360) 2 , (0 oe (/) - 0.0) 2 ] 


< min [(0 O g (f) - 360) 2 , (0 oe (/) - 0.0) 2 ] 

-e- 

o 

cc. 

II 

O 

(33) 

(0o «(/) - 180) 2 

if 0oe («T /) = 360 
if 0oe (to,/) = 180 


CT(/) + 1000C 2 (/) + C 3 (/) 


(34) 


where 


Pk 


10(m/2O) 

0.25 


and Eq. ( 32 ) defines artificial constraint C3 which limits the magnitude of the auto-spectrum calculated 
with the restricted modal model and insures that it is near the measured auto-spectrum. The mode spans 
shown in table 2 on page 12 indicates cutoff frequencies selected. They correspond to those discussed in 
Appendix A. The same table is used for all test conditions. 


IV. Results 


A. Reconstruction of Measured Data 

The first measure of the extent the restricted modal model implemented with numerical optimization pro- 
cedure yields an acceptable set of modal amplitudes is how well the coefficients reconstruct the measured 
cross-spectrum using Eqs.( 22 ), ( 23 ) and ( 24 ). Results are shown in Fig. 6 on page 22. It is clear that the 
procedure discussed is effective in finding a set of coefficients which can construct the measured amplitude 
and phase spectra. 

A second measure of the extent the restricted modal model implemented with numerical optimization 
procedure yields an acceptable set of modal amplitudes is how well the coefficients reconstruct the measured 
auto-spectrum using Eq.( 25 ). Results are shown in Fig. 7 on page 23 The measured auto-spectrunr contains 
noise and should be larger than the reconstructed auto-spectrum. At frequencies less than 180 Hz. where 
only the m = 0 mode exists the results are quite good as expected. At higher frequencies the artificial 
constraint included in the minimization procedure is forcing the reconstructed auto-spectrum to be near the 
measured auto-spectrum. 

B. Model Content 

In order to characterize the model coefficients the modes are normalized by the maximum value. The 
normalized coefficients are shown in Fig. 8 on page 24. 

The mode normalization coefficients are plotted in Fig. 9 on page 25. Fig. 9 on page 25 also shows linear 
curve fits to the normalization coefficients. 

C. Synthesized Cross-Spectra 

The normalized coefficients by mode number for all cases were studied. The uniformity of the plots lead to 
the idea that one could use the mean value of the coefficients to synthesize the cross-spectra at a range of 
angles and shaft rotation speeds. The mean was constructed using the first eight test conditions. Values for 
test conditions at 2304 rpm (N1 CORR) and 2400 rpm ( N1 Corr) were excluded. 

Mean relative coefficients are shown in Fig. 10 on page 26. 

V. Discussion 

A new treatment of the modal decomposition of the pressure field in a PW4098 engine combustor as de- 
termined by two Kulite pressure measurements was developed. A new method using aligned and unaligned 
coherence measurements showed that Kulite data presented is valid. The treatment uses the procedure of an- 
alyzing pressure propagation in the combustor in the same manner that pressure propagation in a compressor 
or fan are analyzed. A similar procedure was used by Karchmer 16 to analyze YF102 combustor coherence 


NASA/TM— 2006-2 14351 


14 


measurements made in a test stand and engine. This procedure for the operating conditions considered uses 
an analysis in terms of waves that propagate in modes such that for waves other than the plane wave the 
mode of propagation is determined by the condition that the frequency of the wave is greater than a cutoff 
frequency. This cutoff frequency is determined by an acoustic eigenvalue equation, wall boundary conditions, 
and a longitudinal boundary conditions that correspond to an infinite duct with no reflections. This type of 
analysis is completely different than one based on the use of duct modes which are solution of the acoustic 
eigenvalue equation, with duct wall boundary conditions and duct exit and entrance boundary conditions. 
The analysis provides evidence the peaks and dips in the measured auto-spectra and cross-spectra magnitudes 
are due to annular duct propagation modes rather than a fluctuating pressure generating mechanism. 

The modal analysis also uses assumption made by Karchmer 16 that equal amplitude circumferential 
clockwise and counter-clockwise pressure waves exist for each circumferential mode, to. However, for modes 
other than the plane wave the new part of the treatment is the assumption that there are distinct frequency 
bands in which the individual modes, including the plane wave mode, overlap such that if circumferential 
mode to and circumferential mode to — 1 are present than circumferential mode to — 2 is not. Consequently, 
in the analysis used herein at frequencies above the first cut-off mode frequency, only pairs of circumferential 
modes are assumed present at each frequency. The measured data up to 800 Hz. was reconstructed. The 
to = 0 plane wave mode was used over the interval [0, 175] Hertz. The m = 0 and to = 1 annular modes were 
used over the interval [188, 328] Hertz. The to = 1 and m = 2 annular modes were used over the interval 
[340, 539] Hertz. The m = 2 and m = 3 annular modes were used over the interval [551, 832] Hertz. 

In agreement with results obtained previously by Karchmer, 16 the results obtained herein indicate that 
for each circumferential mode, the E function for the radial pressure distribution (see Appendix A) 
corresponds to the first root, /i = 0 which has the least radial variation from hub to tip. A new result 
discussed in Appendix A is that the successful use of the same modal span frequencies (see Table 2 and 3) 
over a range of operating conditions for this particular engine (PW4098) suggests that the temperature. T, 
and velocity, v, of the flow at each operating condition are related by c 2 — v 2 = a constant where c is the 
speed of sound, c = jRT/MW . This result is obtained since the mean flow is included in the derivation of 
the annular duct mode equations. 

The comparison of these results for the PW4098 combustor mounted in the engine with similar results 
for the YF102 combustor in an engine and installed in a ducted test rig ( Karchmer 16 )show many similar- 
ities in the pressure field since the pressure field for the PW4098 and the YF102 depends on the first few 
circumferential modes, m = l,m = 2, to = 3... and the corresponding zeroth order radial modes /r = 0. 

Since the YF102 combustor component test rig did not include any turbo machinery this fundamental struc- 
ture of the pressure field seems to be independent of propagating pressure fields from either the compressor 
or turbine or difference tones created by interacting waves from the compressor and turbine. The fact that 
the energy distribution is in frequency bands with well defined spans (see Table 2 and 3) related to particular 
cut-off modes of propagation explains the occurrence of well defined bands in Fig. 4 where the aligned and 
unaligned coherence functions are presented. 

The frequency span of the modes for this case does not resemble that in the YF102 since the PW4098 
larger than the YF102 ( Karchmer 16 ). The energy distribution in the plane wave mode and the first few 
circumferential modes is identified over a range of operating conditions. Due to the difference in size the 
energy distribution in the modes for the PW4098 is again different since the PW4098 is so much larger than 
the YF102 ( Karchmer 16 ). The model has two drawbacks. The first is the assumption that equal amplitude 
circumferential clockwise and counter-clockwise pressure waves exist for each circumferential mode. The 
phase angles shown in Fig 5 would be only 0°,360°, or 180° if this was true. This assumption seems more 
valid for the YF102 data ( Karchmer ( 16 ). The second drawback is shown in Fig. 6 where the reconstructed 
auto-spectrum is far below the measured auto spectrum below 200 Hz where the plane wave propagates. 
Above 200 Hz the difference between the measured and reconstructed auto spectrum is not very large. 

VI. Conclusions 

The treatment of the modal decomposition of the pressure field in a combustor as determined by two 
Kulite pressure measurements developed herein is in reasonable agreement with the measurements. It was 
applied to data from a Pratt & Whitney PW4098 engine combustor over a range of engine operating condi- 
tions. The method works well for the plane wave mode (below 200 Hz). The method gives meaningful results 
at higher frequencies. The method does show it is feasible to use the available measurements to produce a 
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restricted modal model. 


A. Pressure Field in Annular Duct with Flow 


The pressure field, p = p(r, <fi, z, t) of an ideal gas flowing in an annular duct formed by two concentric 
cylinders with inner radius, r a and outer radius, Tb, is a solution of the homogeneous wave equation with 
constant mean density and constant speed of sound. 


1 D 2 p 
S 2 Dt 2 


V 2 p' = 0 


(35) 


where c is the mean acoustic speed, the space coordinates are (r, <t>,z), t is time and D/Dt is the substantial 
derivative, D/Dt = + wjj. 

The combustor will be treated as a hard wall annular duct with a mean axial flow. Propagation of 
acoustic modes in an annular duct with no mean flow is discussed by several researchers . 4, 11,37 The mean 
axial velocity is included in a discussion of axial flow acoustic modes propagating in a cylinder by several 
researchers 7 , 12 ~ 14 > 3 8,39 an d j n an a nnular duct by Envia 40 . Since mode propagation in a annular combustion 
duct with mean flow mean flow is has not been discussed a treatment is given herein. 

Equation (35) in polar coordinates in three dimensions is 


dr 2 


ld_ 

r dr 


Id 2 d 2 1 (d 


r 2 d(j> 2 dz 2 c 2 \dt 


- — I —+u— 

dz 


p = 0 


(36) 


This partial differential equation is converted into an ordinary differential equation if p is taken to be the 
sum of an oscillatory and constant pressure terms, 


p = p 0 {x,t)+p c . 


(37) 


The oscillatory pressure p a is separable into the product of a space-dependent and time-dependent function. 


Po(x,t) = P(x)T(t) 


(38) 


where we assume T = T 0 e lult . 

The space-dependent part, P, is the product three single-argument space functions, 


P{x) = R{r)<S>(4>)Z{z). 


(39) 


where the solutions are $(</>) = cos (f3(f>) or $>((/>) = sin(/30) 
and Z{z) = Z 0 e ikn ). 

Using the convention that primes denote differentiation of a single argument function with respect to 
that argument, the wave equation becomes 


P" 

n 


Iff. 
r R 


1 

7*2 (j) 


or 


1 (T 


V Z' 


— — 37 — - + 2 u—— + u — — = 0 


T Z 


Z 


(40) 


RT 

It 


Iff. 
r R 




R 


1 R‘ 


Z" 


R r R r 2 ^ 


Z' 

'z 


+ ( (f - M 2 )— 2 M^r 

1 Z c 


2 Z "\ 

+U VJ 

= 0 

(41) 

Z' u 2 \ 
Z + c 2 ) 

= 0 

(42) 


where M = uc. 

Now we define the wave number, k = lo/c and choose (— ak ) 2 as the constant equal to the .Z-dependent 
group so that from Eq. (42) 

(l^M 2 )^2¥M^ + ^ + (af) 2 = 0 (43) 

Zj c Zj c 
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and 


(44) 


R" 1 R' 1 , , n2 l2 „ 

— 4 77 j ^ ' ( a ^) + ^ — 0 

R r R r z 

We take as a solution to Eq. 43 
Z(z) = Z 0 e ikQz . Then 

^ = o 

(1 — M 2 ) (1 — M 2 ) 

Equation 45 has two solutions 


(45) 


n+ 

n~ 

We define a new eigenvalue, a, as 


M+ [M 2 + a 2 (l — M 2 )] 1/2 
(1 - M 2 ) 

M - [M 2 + a 2 (l — M 2 )] 1/2 
(1 — M 2 ) 


(46) 

(47) 


a 2 = fc 2 (l — a 2 )r 2 


where rb is the outer radius. Then Eq. (44) becomes 


(48) 


R n 


-R' 

r 


p 2 


L ' 6 


i? = 0 


(49) 


This is Bessel’s differential equation with an eigenvalue, a/rb- Solutions are Bessel functions of the first kind, 
Jp(ar/rb), and Bessel functions of the second kind Yp(ar/rb) which are discussed in Abramowitz. 41 
The space dependence of the pressure is P(r,ip,x) = R(r)'b(tp)Z(x) where 


R(r) = C 1 J f 3(ar/r b ) + C 2 Yp{ar/r b ) = C 1 [J/ 3 (ar/r b ) + QpYp(ar/rb)} (50) 

= C\Ep{ar /r b ) (51) 

$(0) = C 3 sin(/?(/>) + C *4 cos(/30) (52) 

Z(z) = C 5 e ikn+ + C 6 e ikn ~ (53) 


The constants are evaluated by 

• Applying boundary conditions. 

• Requiring the pressure to be a continuous unique function of position. 

• Requiring the pressure to be finite everywhere. 

The wall boundary conditions are determined next. In terms of acoustic velocity, u, where 

du _ 

p— = -Vp (54) 

If the normal velocity vanishes at the boundary so does the normal acceleration. Consequently 

(Vp)„| & = 0 (55) 

On the walls the normal velocity is zero. Thus since if u is zero so is u\ n we have R' (r) = 0 on the walls. As 
a consequence on the wall of the inner shell (r = r a ) and the outer shell (r = r b ) a pressure antinode is set. 
A pressure node is a region in which the air pressure is constant. A pressure anti-node (loop) is a region 
where the pressure is not constant. 

Thus 
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R\r) = 0, 

= o, 


(r = r a ,r b ) 
(<fi = 0 , 2n) 


From equation (52) and (57) 

1 


[/ 3 C 3 cos(/3 </) — f3Ci sin(/3/>)] = 0 , (^ = 0 , 4> = 27t) 

/ 

Then C 3 = 0 and (3 = m where m is the angular index and m = 0, 1, 2.... 

From Eqs. ( 51 ) and ( 56 ) 


(56) 

(57) 

(58) 


c A J ' m {ar/r h ) + C 2 ^Y m (ar/r b ) 


= 0 


r = r a ,r b 


Qm = -7f = 


J'm(a r a/n) 

J'mia) 


C 2 _ J m (ar a /r b ) 
Ci Y^(ar a /r b ) 

Jm( a ) 

Y m( a ) 

Y L{a r a/n) 

Y m(a) 


(59) 

(60) 

(61) 


If the duct has no inner boundary, then r a = 0 and the boundary condition requires that C 2 in Eq. 51 be 
zero, since F/?(0) = — 00 . Then Eq. 51 becomes 


R(r) = CiJ m {ar/r b ) 

The boundary wall condition yields J m (a m i) = 0 when r a = 0. 

Letting a = r a /r b , the eigenvalues of the solution are given by the roots of 


Y'jya^o) 


= 0 


r ytz 0 


(62) 


(63) 


Y m( a mi ) 

Using these values of values of can be determined from Eq. 60 and the radial characteristic 

function of r is R(r)/Ci = + Qmih m (o^^) = 

From Eq. 48 the wave number a is 


a = ±- 


lk 2 - 


>h 2 

) 


(64) 


From Eq. (46) and Eq. (47) the axial propagating mode is 


Z(z) = Z a e 


(1-M 2 ) 


^-{Mk± 


( a M I 2 

fc 2 — (1— M 2 ) - 


1/2 


(65) 


( a mi) 2 


If k 2 is less than (1 — M 2 ) , the traveling wave solution is attenuated with distance. From this rela- 

tionship the frequency,/, of a propagating disturbance must be greater than f cu toff where 


f critical 

f cutoff ) 


C (ami) 
2n r b 


= 1/(1 - M 2 )f cr iti ca l 


( 66 ) 

(67) 


The function f C utoff(M) varies slowly with M when M is low. For example, at M = 0.6, fcutoff = 0 . 8 f critical- 
The reduction of the cutoff frequency by mean flow has been shown experimentally by Mason 42 to follow 
Eq. 67 for Mach numbers less than 0.17. 
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Table 3. Cuttoff Frequencies (T = 1619. 37°K , a = 0.66 , r a = 0.558 m , = 0.837 m , c = 806.7 m/sec 


order m 

zero sequence number p 

(alt) 

J cutoff — 2 nr b ’ 

nW 

mfi 

1 

0 

1.207637 

185.245 

-0.408810 

2 

0 

2.412094 

370.003 

-0.348305 

3 

0 

3.610271 

553.797 

-0.314771 

4 

0 

4.799250 

736.181 

-0.270046 

5 

0 

5.976450 

916.757 

-0.223789 


Table 3 shows, the eigenvalues, Qmjj , and 

( a M ) c 

the cutoff frequencies, / cutoff = for the first 

five order modes and /i = 0 using T = 1619. 37°K 
, a = 0.66 , r a = 0.558 m , = 0.837 m ,and 

c = 806.7 m/sec. These parameters produce a mode 
span chart similar the one shown in Table 2 on 
page 12 which is used in the procedure presented 
herein. These were calculated for M = 0.0. Note 
that these same frequencies can appear at other op- 
erating conditions if when the temperature is in- 
creased the flow also increases so that 


1619 
1 - M 2 


(68) 


This relationship is tabulated in table 4. 

In terms of temperature and velocity we have 


Table 4. Operating temperature verses Mach number 
for constant mode cutoff frequencies 


M 

T 

0.00 

1619 

0.05 

1623 

0.10 

1635 

0.15 

1656 

0.20 

1686 

0.25 

1726 

0.30 

1779 

0.35 

1845 

0.40 

1927 

0.45 

2030 

0.50 

2159 


Cl\J\ - (v 2 /c 2 ) 2 

2 2 
c 2 - v 2 

t 2 


Ci\/l - (iq/ci) 2 
c\ - v\ = A 

Vv 2 + A 
jR/MW 


(69) 

(70) 

(71) 


Consequently, at each operating point it appears that the combustor flow velocity and temperature are 
related by Eq. 71. 

Characteristic d-m/!. («m ; i f - ) functions for radial pressure distributions are shown in Fig. 11 on page 27 
for fi = 0 and m = 1,2, 3, 4, and 5 . 
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Figure 6. Measured and reconstructed cross-spectrum. 
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Figure 7. Measured and reconstructed auto-spectrum and cross-spectrum magnitude. 
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Figure 8. Relative amplitudes, A m (/) = 8 m (/). 
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Figure 10. Mean relative amplitudes for m = 0, m = 1, m = 2, and m = 3. 
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Figure 11. Characteristic E function for radial pressure fuunction for fi = 0 and m = 1,2, 3, 4, and 5. 
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